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Abstract 

We consider both closed and open integrable antiferromagnetic chains constructed 
with the S'C/(AA)-invariant R matrix. For the closed chain, we extend the analyses of 
Sutherland and Kulish - Reshetikhin by considering also complex "string" solutions 
of the Bethe Ansatz equations. Such solutions are essential to describe general mul- 
tiparticle excited states. We also explicitly determine the SU (Af) quantum numbers 
of the states. In particular, the model has particle-like excitations in the fundamen- 
tal representations [k] of SU (Af) , with k = 1 , . . . ,Af — 1. We directly compute the 
complete two-particle S matrices for the cases [1] <g> [1] and [1] (g> [Af — 1]. For the 
open chain with diagonal boundary fields, we show that the transfer matrix has the 
symmetry SU(l) x SU(Af — I) x U(l), as well as a new "duality" symmetry which maps 
I <-» Af — I. With the help of these symmetries, we compute by means of the Bethe 
Ansatz for particles of types [1] and [Af — 1] the corresponding boundary S matrices. 



1 Introduction 



Integrable quantum spin chains are exactly solvable quantum mechanical models of N quan- 
tum spins, of which the Heisenberg model solved by Bethe fl|]-|§ is the prototype. In the 
antiferromagnetic regime, such spin chains can be regarded as lattice versions of correspond- 
ing integrable relativistic quantum field theories. For integrable spin chains, quantities of 
physical interest (spectrum, S matrix, etc.) can be calculated exactly by direct means, 
starting from the microscopic Hamiltonian; while for the corresponding field theories, such 
exact information has been primarily obtained by indirect means, such as the "bootstrap" 
approach [|J, || and semiclassical approximations. 

Quantum spin chains have one spatial dimension and come in two topologies: "closed" 
(periodic boundary conditions) and "open" . The latter exhibit a rich variety of boundary 
phenomena, which - for integrable chains - can be investigated exactly. 

There exists a systematic approach for constructing integrable quantum spin chains, 
called the Quantum Inverse Scattering Method. (For a recent review, see 0.) The basic 
building blocks for constructing closed chains are R matrices, which are solutions of the 
Yang-Baxter equation 

R12W Rxs(X + A') R 23 (X') = R 23 (X') R 13 (X + A') i? 12 (A) . (1.1) 

For constructing open chains, one needs also K matrices, which are solutions of the boundary 



Yang-Baxter equation P]-[10 



i?i 2 (Ai - X 2 )K 1 (X 1 )R 21 (X 1 + X 2 )K 2 (X 2 ) = K 2 {X 2 )R 12 {X X + A 2 )iT 1 (A 1 )i? 21 (A 1 - A 2 ) . (1.2) 
In this paper, we focus on integrable quantum spin chains constructed with the R matrix 
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RW = ^—(^ + w) , (1.3) 

A + 1 

where V is the permutation matrix 

Vx ® y = y <g> x (1.4) 

for all vectors x and y in an A/"-dimensional complex vector space CV- This R matrix is 
SU(J\f) invariant; i.e., 

[U®U,R(\)]=Q (1.5) 
for all group elements U in the defining representation of SU(J\f). 
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This paper has two main sections. In Section 2 we consider the closed integrable antifer- 
romagnetic chain with N "spins" (vectors in CV), which has the Hamiltonian 



N-l 



H 



closed 



nn+l 

+ H 



Nl 



;i.6) 



n=l 



where the two-site Hamiltonian Tijk is given by 

i d 



2dX 



A=0 



(v-r 



;i.7) 



The model is SU(Af) invariant, and the space of states is C® . This a generalization of the 



antiferromagnetic Heisenberg model, which corresponds to the case Af = 2. 

This model was first investigated by Sutherland ]12[ and by Kulish and Reshetikhin [[13 



These authors determined the energy eigenstates and eigenvalues in terms of Af — 1 types of 
roots of a system of Bethe Ansatz equations (BAE). Moreover, they found that the ground 
state corresponds to having Af — 1 "filled Fermi seas" ; and they studied particle-like excited 
states Q corresponding to "holes" in these seas. These analyses considered only real roots of 
the BAE. 

We extend these analyses by considering also complex "string" solutions of the BAE. 
Such solutions are essential to describe general multiparticle excited states. Moreover, we 
explicitly determine the SU(Af) quantum numbers of the states. In particular, we show that 
the Bethe Ansatz state consisting of one hole in the k th sea (k — 1 , . . . , Af — 1) and no 
complex strings is the highest weight of the fundamental representation [k], corresponding 
to a Young tableau with a single column of k boxes, as shown in Figure [I]. That is, the 



>k 



Figure 1: Young tableau with a single column of k boxes, corresponding to 
the fundamental representation [k] of SU (Af) 

model has particle-like excitations in the fundamental representations [k] of SU(Af), with 
k = 1 , . . . ,Af — 1. Finally, we directly compute the complete two-particle S matrices for the 
cases [1] ® [1] and [1] (g> [Af — 1] . In the latter case, the singlet state is described by the Bethe 
Ansatz state consisting of two holes (one each in seas 1 and Af — 1) as well as one string of 
1 Such excitations have been called "kinks" or "spinons". We refer to them here simply as "particles". 
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length 2 in each of the M — 1 seas. Our results for the S matrices coincide with those found 
by the bootstrap approach ||, with no additional CDD factors. 

In Section 3, we consider the open integrable chain constructed with the SU (J\f )-mvaiia.nt 
R matrix ( |1.3|) , together with the J\f x M diagonal K matrices given by JTJJ 



K?J\,&) = diag(a*,... ,^,6*,... , 6 T ) , (U 




where 

a~ = — A , b~ = it;- + A , 

a + = + b + = i(£ + -AT)-\, (1.9) 

for any / G {1,... , M — 1}. The Hamiltonian is given by 

^ 1 d tr K+ o (0,£ + )H m 

H ope „ = ^ K nn+1 + —-^(A.e-) A=0 + trif + ((U+) • (1-10) 

The parameters which may be regarded as certain boundary magnetic fields, break 
the SU{H) symmetry down to SU (I) x SU (J\f — l) xU(l). Moreover, we find a new "duality" 
symmetry which maps I <-> M — I. With the help of these residual symmetries of the model, 
we compute for particles of types [1] and [TV — 1] the corresponding boundary S matrices, 
which describe scattering from the ends of the chain. This is the first direct calculation of 
boundary S matrices for a model whose symmetry algebra has rank greater than one. For 
the case M = 2, we recover the results of Refs. [15 and [|T^ . 



2 The closed chain 

In this section, we consider the closed integrable chain constructed with the 5'[/(A/')-invariant 
R matrix ( |1.3|) . The transfer matrix t(A) is given by 



t(A) =tr T (A), (2.1) 
where To(A) is the monodromy matrix 

T (X) = R 0N (X)---R i(X). (2.2) 
(As is customary, we suppress the quantum-space subscripts 1 , . . . , N of T (A).) The Yang- 



Baxter equation ( |1 . 1|) guarantees the commutativity of the transfer matrix 

[i(A),t(A')]=0. (2.3) 
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The Hamiltonian (|1.6|) is proportional to the logarithmic derivative of the transfer matrix at 
A = 

i d 

\nut\A) 

A=0 

The "momentum" operator P is defined by 



'Hdosed = --^-logt(A) 
2 dX 



(2.4) 



P = -ilogt(O) 
since t(0) is the one-site shift operator. 



(2.5) 



2.1 SU(J\f) generators 

In the defining representation of SU (TV) , we identify the raising and lowering operators 



3 +{k) = e Kk+l , j( fc )=e fc+1)fc , k = 1,... ,/V- 1, 



and the Cartan generators 



Cfc.A: — efc+i^+i , — l,...,jV — 1, 



(2.6) 



(2.7) 



where e k j are elementary TV x TV matrices with matrix elements {ek,i) ab = Sk,a$i,b- These 
generators obey the commutation relations 



(25 k ,i - 5 k ,i+i - 4+1,0 i 



-(0 



{21 



We denote by ■> s n^ the generators at site n, e.g. 



J. 



i k) = l ® . . . ® l® s (fe) ®i 



n = l,... ,iV, 



(2.9) 



and we denote by J 1 * 1 ^, S 1 ^ the corresponding "total" generators acting on the full space 
of states 

N N 



j±(fc) = ^-±(ft) ^ ^(k) _ ^ 



k = 1 ,. . . ,7V- 1. 



?1=1 



n=l 



The SU{J\f) invariance ( |1.5| ) of the R matrix implies that 



t(X) , J ±{k) ] = \t{\) , s {k) 



o 



k = 1,... ,7V- 1. 



(2.10) 



(2.11) 



For future reference, we now relate these generators to the standard Cartan- Weyl basis. 
We set 



E n 



2 

Af-l 



Hi = £ dS® , * = 1,...,7V-1 



(2.12) 
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(with the coefficients /4 still to be determined), and we demand the commutation relations 

M-l 

[Hi , E a ] = aiE a , [E a , £L a ] = £ . (2.13) 

8=1 

The vectors a 1 = (a\ , . . . , a\f_i) are the simple roots normalized to unity a 1 ■ a 1 = 1. One 
finds 

I 



a 



(0,...,0,-^,Ji±i,0,...,0). (2.14) 



From the second relation in Eq. (|2.13| ), we obtain 

1 



a- 



'•li h = 7i 8 i j k . (2.15) 



This implies the important result that \i k = , . . . , , k — 1 , . . . , H — 1 are the 

fundamental weights of SU{H) (see, e.g., |T] 



2.2 Bethe Ansatz and string hypothesis 

We now review the exact Bethe Ansatz solution of the closed SU (J\f) -invariant chain, and 
we use the string hypothesis to recast the BAE into a form which is particularly suitable for 
studying multiparticle excitations. 

Since the operators H c i ose d , P , S^ k > mutually commute, there exist simultaneous eigen- 
states \E , P , S^ k >). The so-called Bethe Ansatz states are the subset of these states which 
are highest weights of SU(M), i.e., 

J +(fc) | ) = 0, k = l,... ,AT -1. (2.16) 



(See, e.g., Refs. M, |18[, |TJ[.) These states have been determined by both the coordinate 



|T2|j and algebraic [13], p0| Bethe Ansatz methods. In the latter approach, the Bethe 
Ansatz states are constructed using certain creation operators (elements of the monodromy 
matrix) depending on solutions {A^} of the BAE 

1 = n e^XP - \$-») n*2(A?>-A< fc >) n e^-X^), 

13=1 (3=1 (3=1 

a = 1,... ,M (fc) , k = \,...,M-l, (2.17) 



2 The remaining states are obtained by acting on the Bethe Ansatz states with the lowering operators 

j-w. 
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where 



A 4- — 

(2.18) 

2 

and = N , = , A^ = A^ = . The corresponding eigenvalues are given by 

E = ~2 S .(1)2 , ' 
a=l Aa + j 

P = -i£lo gei (A«), 

<*=1 

5(fc) = M (H + M (w)_ 2M W i (2.19) 

We adopt the "string hypothesis", which states that in the thermodynamic (N — > oo) 
limit, all the solutions {A^ , . . . , A^ (fc) } are collections of M^ n ^ strings of length n of the 
form (for > 0) 

X (n M = X (n,k) +i ^+l_^ (22Q) 

where j = 1 , . . . , n; a = 1 , . . . , M^ n,fc ^ k — 1 , . . . , A/"— 1; n = 1, , . . . , oo; and the "centers" 
\(ji,k) are rea j total number of A's of type /c is given by 

oo 

M {k) = Y, nM {n ' k) , fc = 1,... ,M- 1. (2.21) 

n=l 

Implementing this hypothesis in the BAE and then forming the product IT/=i over the 
imaginary parts of the strings (see, e.g., 0), we obtain a set of equations for the centers 
\{n,k) gj ven ( U p an overa n sign) by 

n n F nm(^' k) - ^t' k ~ l) ) n n *wa^ - 

m=l /3=1 I lm=l /3=1 



X 



II II ^(Ai"> fc) - A{T' fc+1) ) } , a = 1 , . . . , M< n - fc > , k = 1 , . . . , A/" - 1(2.22) 

m=l p=l 



where 

,2 /\\ „2 



E nm (X) — e| n _ m |(A) e| n _ m | +2 (A) ••■ e n+m _ 2 (A) e n+m (A) , 
Pnm(A) = e_(| n _ m | + i)(A) e_(| n _ m | +3 )(A) •■■ e_( n+m _ 3 )(A) e_( n+m _i)(A) , (2.23) 

and M(">°) = N5 nl , M^) = , A( n > ) = A^ = . 
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Since Eqs. (|2.22|) involve only products of phases, it is useful to take the logarithm. We 
obtain the following important equations for 



h {n > k \\^ k) ) = J^ k) , a = 1 , . . . , M (n ' fc) 

k = 1,... ,J\f- 1 
The so-called counting function h^ n ' k '(X) is defined by 



n 



1, ,oo. 



(2.24) 



i oo Af(™.>fc-i) 



(m,k— 1) n, 

/9 V 



771=1 /3 = 1 
X3 M( m ' fc+1 ) 



oo Mi m ' k ) 

E E '~ Jn " , 

m=l 0=1 



(A - A{T' fe) ) 



where 



—rim (A) 
^nm(A) 



E E ^(A-A^ fe+1) )}, 

m=l /3=1 



,1 — ^n,m)9|n-m|(A) + 2g| n _ m | +2 (A) + • • • + 2g ra+m _2(A) + g n+m (A) 



(2.25) 



- <Z| n -m|+i(A) + g| n _ m | +3 (A) H h g n+m _ 3 (A) + (A) 

and <7 n (A) is the odd monotonic-increasing function defined (for n > 0) by 

q n (X) = 7T + iloge„(A) , -7r < g n (A) < vr . 
Moreover, {J£ l,k '} are integers or half-odd integers which satisfy 

_ j( n , k ) < j(n,k) <- T(n,fc) 

where is given by 



r(n,fc) = i 
J max 2 



M (n,fc) _ i + jr min(m, n) Im^'^ + M (m ' fc+1) - 2M ( 



m,k) 



m=l 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



In deriving the last equation, we assume the prescription || that J^ ,k > — > J^jf' + n f° r 

that {Ji n ' k) } can be regarded as "quantum numbers" of the 
model: for every set {Jj^' k ^} in the range ( |2.28| ) (no two of which are identical), there is a 
unique solution {X^ } ( no two of which are identical) of Q2.24 ). 

Using the string hypothesis, the expressions ( |2.19j ) for the eigenvalues become 



E 



oo M^"- 1 ) 

-EE 

n=l a=l 

oo A/t™' 1 ) 

-E E [^(Ai n,1) ) 

n=l a=l 



71 



S (k) 



J2 n [M (n ' fc_1) + M {n ' k+1) - 2M {n ' k) 
7 



n=l 



(2.30) 



As already noted, = N5 nA , = 0. 

By invoking the string hypothesis, we have transformed the problem of finding complex 
solutions of the BAE fl2.17| ) to the simpler problem of finding real solutions of the equations 
( p.24j ) . We now proceed to discuss the ground state and excitations. 

2.3 Ground state and excitations 

One can argue[] that the ground state is described by only real roots (i.e., strings of length 
1) and no holes. That is, = for n > 1 and = 2Jgg + 1. Hence, = 

N(N — j)/N. Evidently, the ground state lies in the sector where N/N is an integer. 
Moreover, this state has all = 0, and therefore is a singlet of SU(N), as expected for 
an antiferromagnet. Since there are no holes, this state corresponds to a set of N — 1 filled 
Fermi seas. 

Excited states are described by root distributions with holes and (optionally) complex 
strings (i.e., strings of length greater than 1). We let v^' denote the number of holes in the 

j th sea, 

^ = (2j£g + l) - M^) . (2.31) 



Case a: no complex strings 



We first consider the simpler case of excited states with holes but no complex strings. We 
refer to this a. For this case, we obtain from Eqs. (|2.29| ) and (|2.30|) the remarkably 

simple relation 



,C?) 



It follows from Eq. ( [2.12 ) that the Cartan generators have the eigenvalues 

Af-l 



(2.32) 



(2.33) 



i=i 



We conclude that the Bethe Ansatz state with v^> holes in the j th sea and no complex strings 
is a highest-weight state with highest weight fi given by 



3=1 



(2.34) 



3 One argument is based on the observation that a system at temperature T goes to its ground state as 
T — > 0. The thermodynamics for the case Af = 3 was formulated, following [El) - [p3[, in Ref. pj. 
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where are the fundamental weights of SU(N) (see Eq. ( 2.15Q ). In the corresponding 
Young tableau (see Figure [|), the number of boxes in the i th row is equal to Y^f=i" L/ ^- 



Figure 2: Young tableau corresponding to a general irreducible represen- 
tation of SU(M) 

Equivalently, the representation can be denoted by the number of boxes in each column of 
the Young tableau: 



Af-1,... ,J\f - l,J\f -2,. . . ,J\f-2 



V (M-1) v {M-2) „(1) 

In particular, the state with a single hole in the k th sea (i.e., = ^j.) [ s the highest weight 
of the fundamental representation [k], corresponding to the Young tableau shown in Figure 

i 

We label the holes in the range ( |2.28 ) by { } , a = 1 , . . . , v^'. The corresponding 
hole rapidities {X^'} are defined by 

^)(Ag)) = Jt j \ a = l,...,v®, (2.35) 

where /i^ lj ^(A) is the counting function given in Eq. ( |2.25|) with n — 1. 

For N — > oo and for each value of j, the roots {A^ 1 '-"} become dense on the real line, and 
are described by the corresponding densities a^(X) given by 



a 



1 d 



NdX 

Approximating the sums in h^'(X) by integrals using ^ 

Af(i,j) 

N 



1 M (hj) ,00 1 " U) „ 

T7 E 9(X^) « / 9 (\') a^(X') d\'--J2 g{*V>) 



(2.36) 



(2.37) 



leads to a system of linear integral equations 



Af-l 



Af-1 i/( m ) 



E ((* + * ° {m) ) W = + ^ E E W - ^ m) ) 

m=l m=l ct=l 



(2.38) 



4 Here 5(A) is an arbitrary function of A which goes to for A — > ±00. 
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where 



K(A)j m — a>2(X)5 m j — ai(A)(5 mj _i + 5 mj+ i) , 

a„(A) = _1^A) = ^^L^, (2 .39) 
v ; 2tt d\ 2tt A 2 + £ 

as usual * denotes the convolution 

/oo 

/(A - A') g(X') dX' , (2.40) 

and = i/W = o. 

This system of equations is solved by Fourier transforms, for which we use the following 
conventions: 

/(w) = / e^ A /(A) dA , /(A) = — / e-*- A /(w) dw , (2.41) 
and therefore 

a„(cu) = e- n|w|/2 , n>0. (2.42) 

The resolvent 

^w(A) = (5(A) + /C(A)) m 1 m , (2.43) 



has the Fourier transform 12 



£ , \ _ e' H/2 smh (m<H/2) sinh ((A/" - m>)M/2) 
Wl J sinh (JV|u;|/2) sinh (\u\/2) ' 1 J 

where m > = max(m , m') and m < = min(m,m'). The densities cr^(A) are therefore given 

by 



^■)(A) = S W(A) + i £ " Ai m) )^, m - ^ m (A - Ai m) )] , (2.45) 



m=l o=l 

where the ground state densities s^(A) have the Fourier transforms 

^ A _sinh((A/--j)M/2) 



v ; sinh (W> |/2) v ; 



Heuristically, the 1/N terms in Eq. ( J2.45 ) describe the "polarization" of the Fermi seas due 
to the presence of holes. 
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We remark that Eq. ( p. 31 ) can be solved for the integers M^' in terms of the number 
of holes: 



JV-1 



= £7^(0) {-v {k) +N5 Kl 



k=l 



(2.47) 



In particular, the state with a single hole in the k th sea lies in the sector where (N — k) /M is 
an integer. This is a generalization of the fact [Q, [16] that for the Heisenberg chain (A/ - = 2), 
the state with one hole lies in the sector N = odd. 



The energy and momentum are given by 

E : 



M-l vti) 
j=l a =l 



j=l a =l 

where the ground state energy and momentum per site are given by 



(2.48) 



P = 1T 



where ip(z) = -fflogT(z), and p^\X) satisfies 



dX 



J9 (J) (A) = 2ns {j \X), p (j) (0) = -7T 



M-3 



(2.49) 



(2.50) 



In particular, a single hole in the k th sea (which, as we have seen above, is in the fundamental 
representation [k] of SU(J\f)) with rapidity A^ is a particle-like excitation with energy 
ns^(X^) and momentum p( fc )( A 



Case b: including strings of length 2 

In section 2.4, we compute the full two-particle scattering matrix for the cases [1] © [1] and 
[1] © [A/ - — 1]. Each of these tensor products decomposes into a direct sum of two irreducible 
representations [1] © [1] = [1,1]© [2] and [1] © [J\f - 1] = [J\f - 1 , 1] © [A/] , in accordance 



n®n=m © 



Figure 3: Young tableaux corresponding to [1] © [1] = [1,1]© [2] 
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Figure 4: Young tableaux corresponding to [1] © [J\f — 1] = [J\f— 1,1]© [A/] 



with the Young tableaux in Figures and f|, respectively. The calculation of the S matrix 
eigenvalues requires the densities of the two-particle states corresponding to the irreducible 
representations. The two-particle Bethe Ansatz states which are highest weights of [1 , 1] and 
[A/" — 1 , 1] are of the form described above, with only real roots and two holes. Indeed, these 
states have v^) = 2Sj t i and v^> = 5j y \ +5j > j^^i, respectively. However, the two-particle states 
which are highest weights of [2] and [A/] are not of the above form. In order to describe these 
states, we need in addition to two holes also strings of length 2. 

Let us consider the more general case of Bethe Ansatz states with arbitrary values of 
and M( 2j \ with M^ n ^> = for n > 2. We refer to this as case b. For this case, we obtain 
the following generalization of our previous result (|2.32|) for the eigenvalues of the Cartan 
generators: 

S-W = „(J) + M^- 1 ) + M^' +1 ) - . (2.51) 

Moreover, in the expression for J^&2 we can eliminate the dependence on {M^ 1 '^}: 

J%£ = \ [v® + M^'- 1 ) + M^' +1 ) - - l] . (2.52) 



The densities are now given by 

C£ = l 



( i\\)=^\X)-- E ai(A-A^), j = l,...,AT-l, (2.53) 



where er^(A) is the density given in Eq. ( |2.45|) for the corresponding state with the same 
holes but no complex strings. In order to determine the centers A^ 2 '- 7 -* of the 2-strings in terms 
of the hole positions, we recall that h^ 2 '^(X^'^) = J^' • Passing from sums to integrals, we 
obtain the relations 

2*J™ = E E [^-r^H^-Aft^-i + W 

m=l /3=1 

+ E9i(Ai 2j) -A^), « = 1,...,M^), j = l,...,AT-l. (2.54) 

/3=1 
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The energy and momentum are given by the same expressions in Eq. Q2.48 ). 



We now specialize to the cases of interest. We see from Eq. ( |2.51| ) that the two-particle 
Bethe Ansatz state which is the highest weight of [2] has two holes in sea 1 and one 2-string 
in sea 1 (i.e., = 26 jA and M (2,i) = 5 jtl ). Eq. ( ^52]) implies that jg£ = and hence 
J^ 2 ' 1 ' = 0. From the first relation in Eq. Q2.54Q we conclude that the center Ai ^ of the 
2-string is midway between the two holes 



\r> = - (\\ > + , (2.55) 

independently of the value of M. 

Finally, we consider the singlet [A/] two-particle Bethe Ansatz state. This is the state 
with one hole in sea 1, one hole in sea M — 1, and one 2-string in each of the M — 1 seas (i.e., 
!/(?') = S jtX + Sj,M-i and = 1 for j = 1 , . . . , H - 1). We observe that each j[ 2,j) = 0. 



Remarkably, the relations (|2.54j) lead to a linear system of equations for A^ 2 '^ whose resolvent 



is TZjk(0). We obtain for the centers of the 2-strings the result 

AS 2J) = Aj 1} + (Ar *> - A?>) , j = 1 , • • • , Af - 1 • (2.56) 
This result is represented schematically (for Ai > A^ ^) in Figure |[ 

sea X - 1 o >< 



sea 3 x 

sea 2 x 

sea 1 x-e- 



Figure 5: The singlet [A/] two-particle Bethe Ansatz state. Horizontal lines 
represent the M — 1 root distributions forming the Fermi seas, 
circles denote holes in these seas, and crosses mark the centers 
of the 2-strings. 
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2.4 Bulk S matrix 



Following Refs. [25|, [26|, we define the two-particle S matrix St?1®[ fe ] for particles of type [j] 
and [k] by the momentum quantization condition 



f e i P W(.\W)N S \j]®[k] _ ^ |A (i) , A (fc) ) = 0. 



(2.57) 



where the single-particle momentum p^(X) is given by Eq. (p3D|) , and A (i) , A (fe) are the 
corresponding hole rapidities. 

We focus our attention on the cases [1] ® [1] and [1] <8> [A/" — 1], for which the S matrices 
act in CV <8> CV and CV (g> CV, respectively. As already noted, for these cases, the tensor 
product [j] <S> [k] decomposes into a direct sum of precisely two irreducible representations. 
(See Figures ||] and f|.) The Bethe Ansatz states which are highest weights of these irreducible 
representations belong to the cases we denoted a and b, with densities given by Eqs. Q2.45[ ) 
and (|2.53| ), respectively. Specifically, the state corresponding to the completely antisymmet- 
ric Young tableau (i.e., with only 1 column) belongs to case b, and the other state belongs 
to case a. 

We now compute the eigenvalues of Sw®w. Let S a and Sy, be the eigenvalues of 
corresponding to states belonging to cases a and b, respectively. The identity 



(2.58) 



can easily be obtained from Eqs. Q2.36Q and ( p. 50 ). Integrating from — oo to A^ and 
exponentiating, we obtain the relation 



o'P 



U) 



(2.59) 



where h^' j \X^) = J (lj) . Comparing with Eq. (|2~57D , we see that (up to a rapidity- 
independent phase factor) 



S a ~ exp i z2vr/A / (^(A) - s (j) (A)) dX 



(2.60) 



where cx^(A) is given by Eq. (gg5j ) with z/ w ) = 5 m>j + 8 m>k . The integral can be explicitly 
performed using the Fourier-space expression (|2.44| ), as well as the identity |27 







0X ) (I - e~^ x ) 



-fix 







1 - e- x 



dx_ r(/j)r(/j + (3 + 7 ) 

x - ° g T(^ + P)T{jm + 7) ' 



(2.61) 



5 The case [AT— 1] ® [TV— 1] is equivalent to the case [1] ® [1]. Other cases can presumably be treated along 
similar lines; however, they involve higher-dimensional representations of SU(J\f), and the corresponding S 
matrices have more than two distinct eigenvalues. 
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provided Re \i > 0, Re ji > —Re /3, Re // > — Re 7, and Re /1 > —Re ((3 + 7). One finds 




S a = TT ^ ^ ,^ , A = AW-A«. (2.62) 

11 F 1 2l-j-k 1 iAA r 2l+k-j i\\ ' v 7 

In particular, for the states [1 , 1] and [M —1,1], 

^ = JH!111», ^Af»-A«, (2,3) 



r(i + i(-i + zA))r(i- 

r (|-g) r (l + ^(- 1 + ^ ) 
r(l + #)r(i-&(i + iA) 



= ^ . -t^ » A = A?> - . (2.64) 



Although 5 a has been determined only up to a rapidity-independent phase factor, the 
ratio Sb/ S a can be computed exactly: 

S h i2nN (a^ (X)-a^ [X^dX i27r(^ 1 ' j) (-oo)-hi 1 ' 3) (-oo)) _ i2 7r (j^ 1 j) - J^ 1 j) ) 

= II e^A^-A^), (2.65) 

<x=\ 



where a^j (A) is given by Eq. (|2.53| ). In particular, for the states [2] and [A/], we find 



M^r-^) > ( 2 - 66 ) 



^[1,1] ^2 

S W -e 1 fJ 7 (A< 1 "-Ar i) )), (2.67) 



S[M-l,l] VA/ - 

where we have used our results for the centers of the 2-strings ( |2.55| ) and ([2.561) , respectively. 

Finally, we cast our results into matrix form. For the case [1] <g> [1], SU(Af) symmetry 
implies that the complete two-particle S matrix is given by 

fiWi = s [1 , 1] (x)l(i + v) + s m (\) 1 -(i-v) 

= S [1 , 1] (X) i ^, \ = \?-~$\ (2.68) 
iX + 1 



where S^i^A) is given by Eq. ( |2.63| ). Moreover, for the case [1] ® [Af — 1], 



s imw-i] = s w ^ ltl] (x) (1 - p m ) + s w (x) p m 

= V-i,i](A) U-J^Pm), \ = X$ ) -\? f - 1) , (2.69) 
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where Py^ is the projector onto the one- dimensional subspace [TV], namely, Pysi) = j/Pvi\ 
and S , [jv_i i i](A) is given by Eq. (|2.64|) . These results agree with those found using the 
bootstrap approach by Ogievetsky et al. without additional CDD factors. 



3 The open chain 

We turn now to the open integrable chain constructed with the SU (/V) -invariant R matrix 
(fyD and the diagonal K matrices ( |1.8|) . Although this R matrix does not have crossing 
symmetry for M > 2, it does have the property |28|] 

RuiXY 2 )' 1 )' 2 ) oc M 2 R 12 (\ + 2p) M 2 X , (3.1) 



with M = 1 and p = Af/2. One can therefore prove ||, |2{| the commutativity of the 
transfer matrix (A , £ T ) given by [] 

t (0 (A,e T ) =tr K ( | )0 (A,e + ) T (A) ^ ( 7 )0 (A,£-) T (A) , (3.2) 

where To(A) is the monodromy matrix (|2.2j ) and To(A) is given by 

f (\)=R w (\)---R m (\). (3.3) 

The Hamiltonian ( |1.10|) is related to the derivative of the transfer matrix at A = 



1 d 



i d 



Ko P en= AC - ^ %(A ) e T )^ - i ^logtr^ ( t ) (A,e + ) (3-4) 



A=0 



3.1 Symmetries of the transfer matrix 

The SU(N) invar iance ( |1.5[ ) of the R matrix implies that 

U 2 R 12 (X) U\ = U\ R 12 (X) U x (3.5) 

for all U G SU (/V) . The LHS can be regarded as a "quantum-space" transformation, while 
the RHS can be regarded as an "auxiliary-space" transformation. The quantum-space oper- 
ator IA defined by 

U = U l U 2 ---U N (3.6) 



6 We take into account apparent typos in their Eqs. (2.19) and (2.23) 

7 The more general transfer matrix i(; +j ;_)(A , £ T ) constructed with j(A,^) also forms a one- 
parameter commutative family. For simplicity, we consider here the special case l + = l_ = I 
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therefore has the following action on the monodromy matrices 

UT Q {\)tf = U$T (\)U , 

Uf {\)U ] = U$f Q (\) U Q , (3.7) 
and the transfer matrix transforms as follows: 

Ut (t) (\,&)rf = ti {(u K+(\,£ + ) U^) T (X) (UKy(\,Z-) W)ofo(X)} . (3.8) 
For U e SU(l) x SU(N-l) x U(l), evidently 

UKfc(\,t w )Ut = Kfo, (3.9) 

and therefore 

Wt (0 (A ) ^)W t =t( (A,^). (3.10) 

That is, the transfer matrix has the invariance SU(l) x SU(J\f — I) x U(l). In particular, it 
commutes with all the SU(Af) Cartan generators 

t (0 (A,£ T ),S«] =0, fc = l,...,jV-l. (3.11) 

The transfer matrix also has a less evident — but very useful - - "duality" symmetry 
which maps I <-> M — I. This symmetry originates from the simple fact that under the 
transformations 

£- 

-e++AT, (3.12) 
the elements ( |1.9| ) of the matrices transform into each other: a T «-> — 6 T . Therefore, 

^ ( T)(A,^)^-^(A,^), (3.13) 

where 

tfg(A,£ T ) =diag(b^,... ,6 T ,a^,... , a T ) . (3.14) 

l M-l 

Notice that and Kfc are related by a cyclic permutation. Thus, there exist matrices 

U { i) e SU{N) such that 

gffl ^(T)(A,^)^) = ^- ( A '^) = -^ (A,r T ) , (3-16) 



8 Evidently, the N y. N matrices E/(j.fc) ( j , fc = 1 , . . . , A/" with j 7^ fc) corresponding to rotations by 7r/2 
in the (j, k) plane are elements of SU (N) . These matrices have the property of permuting the j th and k th 
elements of an arbitrary diagonal matrix D = diag (di , . . . , d^f): 

U(j,k) diag (efi , . . . ,dj , <4 , . . . , d^) U^ k) = diag (di , . . . ,d k ,... ,dj , dtf) ■ (3.15) 

We choose the matrices Lfy) to be products of matrices of the type Ufj^y Note that Eq. fl3.16| ) does not 
uniquely determine Um ■ 
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where 

e + = -e + +A/-, 

V = M-l. (3.17) 

Correspondingly, we obtain the desired "duality" transformation property of the transfer 
matrix 

%%)(^tKo=%')( A ^)> ( 3 - 18 ) 

where 

= i U(i) 2---U$) n ■ 

Notice that the square of this transformation is the identity. The transfer matrix is "self- 
dual" for £_ = £_ = 0, = £V = AT/2, and Z = Z' = Af/2. 



3.2 Dual pseudovacuum 

The algebraic Bethe Ansatz is usually implemented with the pseudovacuum oom, where 



/ \ 



(3.19) 



V o y 



We shall also make use of Bethe Ansatz states constructed with the "dual" pseudovacuum^ 

J = U(i) = ^ (A T) (3.20) 

in order to compute boundary S matrices. We can obtain the corresponding Bethe Ansatz 
equations (BAE) with the help of the following 

Lemma Let g be a transformation on the parameters I and £ T of the boundary K matrices, 



9 



(3.21) 



9 We use here the matrix Un\ described in the previous footnote. A different choice of can lead to a 
different "dual" pseudovacuum. Nevertheless, the results given in Section 3.4 for the boundary S matrices 
do not depend on this choice. 



18 



which squares to the identity, i.e., g 2 = 1. Furthermore, letUm be a unitary transformation 
on Cff such that 

U {l) t {l) (X,^) W$> = t {n (X,Q (3.22) 

U.(i) uj = uj' . (3.23) 

Then the BAE for the transfer matrix t(z)(A , £ T ) with the pseudovacuum uj' are the same as 
the BAE for the transfer matrix t^)(X,^' T ) with the pseudovacuum uj. 

Before giving our general proof, it is instructive to examine a more explicit proof for the 
special case j\f = 2, which was first considered by Sklyanin ||. For this case 1 = 1' = 1, and 
we therefore suppress the label / and write the transfer matrix as t(A,^ T ). Moreover, we 
take U as in Eq. (|3.6|) , with 

U = | I | ■ (3.24) 



Then condition ( [3.22j ) is satisfied, with £' as in Eq. ( |3.17|) . We consider the pseudovacuum 

uj = uJ(i), and thus uj' = U = ujq). The algebraic Bethe Ansatz leads to the fundamental 
result 

t(X,£ T )\{X a },^) = A{X,{X a },^) |{A Q },U, (3-25) 

where 

| {A Q } , £_) = B(\i , £_) • • • B(X M (3.26) 

and {A Q } are solutions of the BAE with pseudovacuum uj. Multiplying both sides of Eq. 

( gg| ) with U, and using condition ( ggg ) as well as the fact U £(A W f = C(X,£'_), we 
see that 

f(A,£y |{A Q },0 / = A(A,{A Q },e T ) |{AJ,0', (3-27) 

where 

|{A Q } , £_>' = C(A X , £_) • • • C(A M , £-) • (3-28) 
Changing £ T — > £' in Eq. ( |3.27| ), we obtain 

t(\,&) |{A / J,e-)' = A(A,{A / Q },^) \{K>},t-Y, (3-29) 
where {A^} satisfy the same BAE as {A Q }, except with £ T — > C T - 
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We consider now the general case. The nested algebraic Bethe Ansatz leads to the result 



%(A,£ T ) | ) = A (0 (A,£ T ) | ), (3.30) 

where the eigenstate | ) is constructed with the pseudovacuum u>. Multiplying both sides by 
Um and using condition ( |3.22j ) gives 

t (n (\,e T )U {l) |> = A (I) (A,e T ) U {1) | ). (3.31) 

Changing £ T — > £L and I — > I', we obtain 

t {l) (X,^) | )' = A (n (A,^) | >', (3.32) 

where | )' is constructed with the pseudovacuum oj' = Uq) uj. Notice that the eigenvalue of the 
transfer matrix is A^)(A Recalling (see, e.g., 0) that the BAE are precisely the condi- 
tions that the eigenvalues have vanishing residues, we conclude that the BAE corresponding 
to the pseudovacuum to' are the same as the BAE corresponding to the pseudovacuum uj, 
except with £ T — ► £' and I — > I'. This concludes our proof of the Lemma. 



3.3 Bethe Ansatz and multihole states 



The eigenstates of the transfer matrix t(^(A,^ T ) with the pseudovacuum u>^ have been 
constructed in Refs. ]TJ|], [j3"D| . The corresponding Bethe Ansatz equations (BAE) are given 
byQ 

eH-+i(*8>) e- ( 2 f+ -2Ar+/)(AL /) ) 8 l)k + (1 - 6 l)k ) 



M( fe+1 ) 
/3=1 



/9=1 



a = 1 



k = 1 



1. 



(3.33) 



As before, = N , M<M = , A^ 0) = A ^ = . The requirement that solutions of 

the BAE correspond to independent Bethe Ansatz states leads to the restriction A^ > 0. 
For later convenience, we restrict > \{N — 1), 6+ > A/" — \. (See Eqs. (]3.42|) and ( |3.43| ) 
below.) 



10 c 



■'Starting from Eq. (17) in [Hfl , we make a shift of variables /x^ — > (i^ — 4^; and we then take the 

(k) (k) 

isotropic limit by making the redefinitions 7 = irj, £ T — > 177^, = r/X^ , and then letting 77 — > 0. 
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Since we need to consider only one-particle states in order to compute boundary S matri- 



ces 



161 , we restrict our attention here to real solutions of the BAE, i.e., no complex strings. 



In terms of the counting functions hi^(X) defined by 



<) } (A) 



2tt 



{91(A) + [-q 2 z_+i(\) + g 2 £ + -2^+«(A)] 4,/ 



(3=1 



E 

(3=1 



Y: [ 92 (A-A«) + g2 (A + A 



> 



+ E [^A-A^V^A + A^)]} 



(3=1 

the BAE take the form 



k = l t ... M- 1- 



(3.34) 



(3.35) 



Although ultimately we focus on one-hole states, it is convenient to first consider the 
more general case of multihole states. As in the previous section, we let v^> denote the 
number of holes in the j th sea, and we define the hole rapidities {A^} by 



a = 1 v 



U) 



In the thermodynamic limit, the roots are described by densities t7$(A) given by 



The sums in fom (A) can be approximated by integrals using (see, e.g., flT5| ) 



1 M<J) 

jr E 9^) 

a>=l 



1 vK ~ 1 

«7(A') dA' - - E 0(A<?) - 2^0(0) • 



For the symmetric density c^m S (A) defined by 



<$.(A) 



<$(A) A>0 
<$(-A) A<0 



we obtain the system of linear integral equations 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



E ((* + *0 im * 4T } ( A ) = 2 °i(A)^i + ^{« 2 (A) + ax (A) (-1 + S jA + fa-i 



m=l 



+ (-a 2 ^ +i (A) + a 2S+ _^(A)) * ifI + E E fc(A - Ai m) ) jm + /C(A + A^);™) } , 

m=l a=l 

j = l,...,X-l, (3.40) 
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where /C(A) jm is defined in Eq. ( f2.39|) . The solution is given by 



a 



(J) 
(0 ■ 



i JV-1 

2s® (\) + 6<r$(\) + ^{ E (^-im * [o 2 + en (-1 + <5 m)1 + 5 m ^-i)]) (A) 



(A) 



m=l 



+ "E E - Ai m) )5,, m - 7^ m (A - A>*>) + (A<T> 



m=l a=l 

,Cj')/ 



where the quantity Sa^(X) defined by 



-2AT+Z 



(A) 



Ai m) )] } , (3.41) 



(3.42) 



has the dependence on the boundary parameters and 7£j m (A) is the resolvent, which has 
the Fourier transform ( |2.44| ). 

We shall also need the densities a'M\(X) corresponding to the dual pseudovacuum u' 
given by Eq. ( |3.20|) in order to calculate boundary S matrices. According to the Lemma, 
the BAE with the pseudovacuum uj' are given by Eq. (|3.33|), except with £ T — > £L and I — > /'. 



It follows that the corresponding densities cr'^\(X) are given by Eq. (ft.41|) , except with 



^(!) } (A) = (Rjjf-i * (a 2i __ M+ i - a 2 ^ + -M+i)) (A) . 



(3.43) 



As previously noted, the K matrices in the transfer matrix tm(X,^) break the bulk 
SU(J\f) symmetry. Hence, strictly speaking, states should be classified according to the 
unbroken symmetry SU(l)xSU(J\f—l) xU(l). However, we expect that at points of the chain 
that are far from the boundary, the effects of the boundary should be "small" . In particular, 
in the bulk, multiparticle states should "approximately" form irreducible representations of 
SU (A/ - ), as discussed for the closed chain in Section 2. Therefore, we shall continue to classify 
bulk multiparticle states by SU (J\f) quantum numbers. 

Consider now the special case of one-particle states of type [1], which form an M- 
dimensional representation. For the Bethe Ansatz state constructed with the pseudovacuum 
LO{k\ having one hole in sea 1, the Cartan generators have the eigenvalues R 

S® = 5 jtk - 8 j>k - x . (3.44) 

This state is represented by the vector 

/0\ 



k 



th 



(3.45) 



11 In particular, — 6j t i for k = 1, which is consistent with Eq. 
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which is the eigenvector of the matrices {s^'} given in Eq. ( |2.7|) with the eigenvalues ( 3.44Q . 

Similarly, the one-particle states of type [Af — 1] = [1] also form an Af- dimensional 
representation. The Bethe Ansatz state constructed with the pseudovacuum loha having one 
hole in sea Af — 1 has the eigenvalues 



This state is represented by the vector 

/ \ 



(3.46) 



(AT+l-k) 



th 



(3.47) 



v o y 



keeping in mind that the Cartan generators are now represented by the matrices {— s^'*}. 
The highest weight of [1] is the negative of the lowest weight of [1]. 



3.4 Boundary S matrices 

We define the boundary S matrices ^ for a particle of type [j] , in analogy with the bulk 
S matrix, by the quantization condition |fL5| , |T6| 

(e^ U) ^)» S + m S w w - l) |A«> = 0, (3.48) 

where pV'(\) is defined by Eq. ( |2.5(J| ), and is the hole rapidity. There is an identity for 
the open chain which is analogous to the one given in Eq. ( 2.58j) for the closed chain: 

iix pij){x) + ^ (A) - 2s °' )(A) = ivi\ h v (A) (3 - 49) 

For simplicity, we focus our attention on the cases [1] and [Af — 1] = [I] , for which the S 
matrices act in CV and CV, respectively. 

We first treat the case [1]. The SU(l) x SU(Af — I) x U(l) invariance of the transfer 
matrix implies that m are diagonal Af x Af matrices of the form 

Sfo {1] = diag (a^ , . . . , ofa , , . . . , . (3.50) 

^ v ' " v ' 

I Af-l 

Choosing the state [X^) in Eq. ( |3.48| ) to be the Bethe Ansatz state constructed with the 
pseudovacuum cj(i) having one hole in sea 1, we see that (up to a rapidity-independent phase 
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factor) 



a (i) «(/) ~ 



exp j^vriV ) (<7$(A) - 2s (1) (A)) dX 



(3.51) 



where o"»V(A) is given by Eq. fl3.41| ). Using the fact 



Ad) 



7£ (A — X^) + TZ(X + A (1) 



dA 



2ft (2A) rfA 



(3.52) 



and with the help of the identity Q2.61 ) and the duplication formula for the gamma function 

1 



2 2 '- 1 J •(.;•) V(.r + - 



7T2r(2a;) 



(3.53) 



we find 



a 



(0 



S (A< 



i) 



r(> (e- + i-| + a«))r(^(e_+^-|-zA« 



AW))r(^(e-+^-i 



a 



(0 



^o(A (1) ); 



r(i (e + -AT + /-i-zAW))r(^(e + 



zAW 
zA« 



r(£ (£ + -^+M+*A«))r(£(e 

where the prefactor So (A) is given by 

, fl . _ r(^(^-i)-»A))r(^(^ + »A 
U r(i(|(^-i ) + zA))r(i(^-zA 



(3.54) 



(3.55) 



Moreover, choosing the state \X^) in Eq. ( J3.48 ) to be the Bethe Ansatz state constructed 
with the dual pseudovacuum uj' = having one hole in sea 1, we obtain the relation 



"J) a (i) 



exp LirN £ ( ' (*$\\) - a§(X)) dX 



(3.56) 



Note that 



<$>(A)-<$(A) = ^(A)-^«(A) 

(a2e_-l(A) - a 2 £ + -2N+2l-l 



(3.57) 



where 5crLV and Sa'/jl are given by Eqs. ( |3.42| ) and ( |3.43| ), respectively. We conclude 



(0 



-e 2 ^_i(A«) 



^| = -e^-^+a-iCA^) 



(3.58) 
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where we have resolved the sign ambiguity by demanding that the S matrix be proportional 
to the unit matrix for A'- 1 -' = 0. 

Finally, we consider the case [Af — 1]. The boundary S matrices S^ nj-_x] are diagonal 
M x TV matrices of the form 

5 (T) [N-i] = diagfaf^ , . . . , aj t) , 0* , . . . , ^) . (3.59) 

S v ' V v ' 

I N-l 

For this case we must consider one hole in sea TV — 1 . Noting that 

pfo ~ exp U2nN jf «~ 1} (A) - 2 S ^~ 1 )(A)) d\\ , (3.60) 
(see Eq. fl3.47| )), we obtain 

^ = 5 ' Q (AC^-)^ r + Z + - X ) - ^^--))) F ( e _ + |(7V- - 1) + iA^-^))) 



r (£- + z + 5^ " !) + ^A^- 1 ))) T (e_ + |(AT - 1) - iACV-D 

00 r a fe + + i - m+ 1) - iA(^-D^ r a fe + - w+ d + 



- (e + + z - |(jv+ i) - zA^-^jj r ^ (e+ - ic/v 

(3.61) 

where So (A) is given in Eq. ( |3.55| ). Moreover, 

^ = -c ae _ +a ^_ 1 (A^- 1 )) , 

5? = -e 2f+ -Ar-i(A^- 1 )) . (3.62) 

4 Discussion 

We have shown how to describe general multiparticle states of the antiferromagnetic SU(Af) 
chain, in particular their SU(J\f) quantum numbers, within the framework of the Bethe 
Ansatz/string hypothesis. The picture which emerges is a rich generalization of the M = 2 
case P|) HI- The ubiquitous appearance of the kernel (l + K(u)) , which is characterized 



(see, e.g., |[20|| ) by the SU(Af) Dynkin diagram, is noteworthy. Moreover, we have computed 
both bulk and boundary scattering matrices for particles of types [1] and [A/" — 1]. It should 
be possible to extend this analysis to particles of any type [k]. 

We have identified a "duality" symmetry of the open chain transfer matrix with diagonal 
boundary fields, which plays an important role in our computation of boundary S matrices. It 
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may be interesting to investigate further the "self-dual" case. We expect that this symmetry 
is also present for the boundary A^_ l Toda theories with diagonal boundary fields. Whether 
such symmetries persist for nondiagonal boundary fields is also an interesting question. 

The "mixed" boundary condition case / + 7^ Z_ merits further investigation. One expects 
that the boundary S matrix for one end of the chain should be independent of the bound- 
ary conditions at the other end. However, for this case, the unbroken symmetry group of 
the transfer matrix is smaller, and therefore, the arguments presented here require further 
refinement. 

Since the groups SO(2Af) and Ejj-, like SU(J\f), have simply- laced Lie algebras, it should 
be possible to treat the corresponding integrable chains in a similar fashion. Finally, we 
note that the boundary S matrix calculations presented here can be generalized ]31[ to the 
trigonometric case, i.e., the open chain constructed with the A$_ 1 R and K matrices. 
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